Quantum Brownian motion for periodic coupling to an Ohmic bath 
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We show theoretically how the periodic coupling between an engineered reservoir and a quantum 
Brownian particle leads to the formation of a dynamical steady state which is characterized by 
an effective temperature above the temperature of the environment. The average steady state 
energy of the system has a higher value than expected from the environmental properties. The 
system experiences repeatedly a non-Markovian behavior - as a consequence the corresponding 
effective decay for long evolution times is always on average stronger than the Markovian one. We 
also highlight the consequences of the scheme to the Zeno-anti-Zeno crossover which depends, in 
addition to the periodicity r, also on the total evolution time of the system. 

PACS numbers: 03.65 Yz, 03.65 Xp 



I. INTRODUCTION 

Fundamental research on open quantum systems has 
traditionally focused on interactions between the reduced 
system and its natural environment [J 0] ■ Recent the- 
oretical and experimental developments in engineering 
the properties of the environments open a new avenue 
by creating artificial controlled reservoirs with whom the 
reduced system interacts 0, 0, [!]• These developments 
make possible, e^-, to observe directly the appearance 
of decoherence [J, Q and to study non-Markovian dy- 
namics of open systems in new regimes [(i, J,]. Moreover, 
reservoir engineering and environment mediated schemes 
allow the indirect control of quantum systems [1, Q - 
the quantum control being the essential ingredient in the 
development of quantum simulators. Generally speak- 
ing, the importance of the open system studies stems 
from their central role in understanding the quantum- 
classical border (lol . [Til , [l^ and from the role of deco- 
herence as an obstacle for creating quantum information 
processors [H, • Reservoir engineering provides a so 
far largely unexplored way to approach these issues which 
have both fundamental and applicative character. 

We focus on the reservoir engineering schemes for a 
paradigmatic open system model, i.e. quantum Brown- 
ian motion (QBM) [1, 2, 15, ij, [H [11, ilB,] , which has 
a wide variety of applications ranging from quantum op- 
tics [13 and nuclear physics [2l[ to chemistry ^2^. An 
earlier study has revealed various types of non-Markovian 
dynamics of QBM for structured reservoirs |6|, i.e. time- 
independent engineered reservoirs psj . Here, we con- 
centrate on a scheme where the system is forced repeat- 
edly to experience non-Markovian behavior due to the 
structured reservoir. We show that this corresponds 
to the case when the engineered reservoir is periodi- 
cally switched off and on - and can also be seen as an 
initial step towards more sophisticated time-dependent 
reservoir engineering techniques. We call this switch 
off-on scheme as shuttered reservoir and consider the 
cases where the duration of each period is on the non- 
Markovian time-scale [l^l . 

So far time-dependent reservoir engineering schemes 



have been considered mostly for few-level systems [25j. 
Here, we give an example of the effects of periodic system- 
environment interaction for a more complicated quantum 
system, namely harmonic oscillator. Our results can find 
applications, e.g., in the field of linear quantum ampli- 
fiers. In that context, the switching off-on action corre- 
sponds to shuttering the pumping lasers, which may lead 
to the interesting effects such as the enhancement or the 
reduction of the gain of the amplified field ^26] . 

It is worth noting that artificial engineered reservoirs 
has been recently created for trapped ion system by ap- 
plying electric noise to the trap electrodes and by using 
laser light jjj, Thus, one can think of creating shut- 
tered reservoirs by shuttering the noise which is applied 
to trap electrodes or by shuttering the laser light. 

The theoretical description of the system dynamics, 
that we consider, is based on the recursive use of the 
appropriate master equation and allows analytical solu- 
tions to be found. For a small number of shuttering peri- 
ods, one can control the appearance of the quantum Zeno 
(QZE) and anti-Zeno effects (AZE) [H, ^ US] for the 
reduced system by tailoring the properties of the envi- 
ronment |7|]. In this paper we concentrate on the inter- 
mediate time dynamics (intermediate number of periods) 
and on the steady state properties of the reduced system 
for long time (large number of periods). We show that 
for a large number of cycles, the system always reaches a 
steady state which is not in thermal equilibrium with the 
environment at temperature T but can be described by 
an effective temperature T^ff > T and by the correspond- 
ing thermal state. The steady state value of the average 
system energy is given by the ratio of the effective time 
averaged diffusion and dissipation constants. 

The scheme has also interesting fundamental conse- 
quences for the appearance of QZE and AZE. As was 
shown in Ref. 0], for a small number of cycles, one can 
control the appearance of the QZE or AZE by chang- 
ing the environment parameters. In this paper we show 
that for large number of cycles the system always expe- 
riences an AZE. In other words, the decay of the system 
is stronger than the Markovian one for large number of 
periods. Thus, for short evolution times the system may 
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experience QZE which for long times turns to AZE. This 
demonstrates that the appearance of quantum Zeno or 
anti-Zeno effects in QBM depends also on the total evo- 
lution time of the system (or number of measurements) 
and demonstrates the richness of Zeno dynamics in QBM 
when compared to more simple systems, e.g, two-level 
atoms. 

The paper is organized in the following way. Section 
im introduces the quantum Brownian motion model for 
a harmonically bound particle and introduces the ba- 
sic properties and parameters for the structured reser- 
voirs we use. These form the framework for Sec. IIIII in 
which the formal treatment of open system with shut- 
tered reservoir is carried out leading to the concept of 
recursive master equation. The results for the heating 
function dynamics are presented in Sec. IIVI and the dis- 
cussion concludes the paper in Sec. |Vl 



II. QUANTUM BROWNIAN MOTION 

We consider a harmonic oscillator linearly coupled with 
a reservoir modelled as an infinite set of non-interacting 
oscillators i, Hi [M [13, E [ll- The dynamics of 
a damped harmonic oscillator is described, in the secu- 
lar approximation, by means of the following generalized 
master equation in the interaction picture [y, |3l| 



dp{t) ^ m+lii) 
dt 2 



\2a) p{t)a — aa) p{t) — p{t)aa) 



(1) 



In this equation, a and are the annihilation and cre- 
ation operators, and p[t) is the reduced density matrix of 
the system harmonic oscillator. The non-Mar kovianity is 
characterized by the time dependence of the coefficients 
A(t) and ^{t) appearing in the master equation - these 
are known as diffusion and dissipation coefhcients, re- 
spectively [gI. [3l|. 

The diffusion coefficient appearing in the master equa- 
tion ([T|), to second order in the dimensionless coupling 
constant g, can be written in the units of h as 
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(1/r) sin(c^oi)]}, 
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(2) 



where the assumption of the high-temperature reservoir, 
n = k-QT/uJo :s> 1 has been used. The dissipation coeffi- 
cient "f{t) can be written as 
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7ft)= 2,1 



-i{ujot)+rs\a{ujQt)]] . (3) 



Above, r = uid^o is the ratio between the environment 
cut-off frequency Wc and the oscillator frequency wq, 



the Boltzmann constant, and T the temperature of the 
environment. 

A commonly used environment for open quantum sys- 
tems is described by an Ohmic reservoir spectral density 
with the Lorentz-Drude cut-off [1] 

•^(^) = (4) 

The spectral distribution is given by 



I{uj) = J (uj)[n^(uj) + 1/2] 



(5) 



where ne{uj) is the occupation of the environment mode 
at frequency uj and Eq. ^ has been used. For high T, 
Eq. ^ becomes 



2kBT 



(6) 



The central parameter r = oJc/ojo describes how on- 
resonant the oscillator is with the reservoir. When r > 
1, intensive part of the environment spectrum overlaps 
with the oscillator frequency and the decay coefficients 
A{t) ± '-f{t) > for all times. Consequently, the mas- 
ter equation is of Lindblad-type [H, |33]. When r < 1, 
the most intense part of the environment spectrum lies 
in the small frequency range and the on-resonant inten- 
sity is small. The decay coefficients A{t) ± 'y{t) acquire 
temporarily negative values and the master equation is 
of non-Lindblad-type [27j . 

The full solution of the master equation ([T|) can be 
found e.g. in Refs. [1, [3l|. In what follows we study 
the time evolution of the heating function {n{t)) with 
n ~ a^a quantum number operator. The dynamics of 
{n{t)) depends only on the diffusion coefficient A{t) and 
on the classical damping coefficient 7(t) [3l[. Further- 
more, the quantum number operator n belongs to a class 
of observables not influenced by the secular approxima- 
tion [ll,[33. 

The solution for the heating function, valid for all times 
and all initial states, is given by 0, [23| 

{nit)} = e-i^(*)(n(0)) + ^ (e-^^W - l) + ArW, (7) 

where the quantities Ar(i) and T(t) are defined in terms 
of the diffusion and dissipation coefficients A(t) and 7(t) 
respectively as follows 



T{t) = 2 f -fih)dti 
Jo 



(8) 



Ar{t) = e-r(*) / e^(*i)A(ti)dti. (9) 
Jo 

For short times the system experiences non-Markovian 
dynamics and for long times, after Markovian heating, 
the system reaches a thermal steady state with its envi- 
ronment. This typical dynamics is modified in the con- 
text of the scheme presented in the next section. 
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III. SHUTTERED RESERVOIR AND 
RECURSIVE MASTER EQUATION 

We are interested in the reduced system dynamics 
when its engineered artificial reservoir is switched off and 
on in a repeated manner with period r - for this pur- 
pose we use the term shuttered reservoir. At the end 
of each period the engineered reservoir is switched off, 
and switching on initiates a new time evolution period 
for the system. Consequently, each switching off and on 
process will reset the correlations between the system and 
the environment. For an experimental realization of the 
scheme, one can, e.g., uncouple the system from the en- 
gineered bath in a periodic way. The duration of the free 
evolution period has to be chosen long enough to avoid 
correlations between different system-bath coupling peri- 
ods, i.e., the duration of the free evolution has to exceed 
the bath correlation time. 

By resetting the correlations between the system and 
the environment, the system oscillator is forced to expe- 
rience non-Markovian dynamics in a periodic way. Con- 
sequently, the system remains repeatedly in the non- 
Markovian regime. Following this scheme, the solution 
for the reduced density matrix and the heating function 
dynamics can be obtained by the recursive use of the 
master equation given by Eq. ([T]). In other words, the 
master equation is solved for each evolution period by 
using as an initial state of the system the one obtained at 
the end of the previous period. In the case we consider, 
the free evolution between system-reservoir interaction 
periods can be neglected in the calculations (33 |. 

This recursive use of the master equation ^ leads, 
assuming an initial Fock-state, to the following equation 
for the density matrix of the reduced system: 



dp{t) 
dt 



= 71 W 
+ 7-1 M 



ap{t)a) 



p{t)a — —aa^p{t) — —p{t)aa' 



(10) 



Here, time t = mr where m indicates the number of the 
period. The decay coefficients are written as 

j^{r)^- r dt'[A{t')+^{t% (11) 



7-i(r) = - f dt'm')^^[t')]. (12) 

T In 



and 



It can be shown that the master equation (fTO|l corre- 
sponds formally to the case where the reduced system 
dynamics is modified by periodic non-selective measure- 
ments. Thus, our shuttered reservoir and initial state 
are equivalent to performing periodically non-selective 
measurements of the energy of the oscillator [s^. For 



a formal study of non-selective measurements and mas- 
ter equations, see Ref. (36j . 

We are mostly interested in the case where the oscil- 
lator is initially in its ground Fock-state \n — 0) (such 
as laser cooled single trapped ion) though the results for 
the steady state presented in the following are valid for 
all initial states. For an initial ground state, the solution 
of the master equation (jlOp . in terms of the quantum 
characteristic function (QCF) x reads 



XtiO = exp 



(nit)) 



(13) 



where the heating function is given by 



{n{t — mr)) — 



l-e-nr) 2" 



(14) 



The solution given in Eq. is obtained by using cer- 
tain algebraic properties of the superoperators of the cor- 
responding generalized master equation, for more details 
see Refs. @,|2l|3l|. Moreover, the QCF of Eq. ([13]) cor- 
responds to a thermal state at all times (c.f. Ref. [37j) 
and its time-dependence is given by {n{t)). 

Eq. p4)) reveals that there exists a steady state value 
for the heating function, (n)s, and its value can be easily 
obtained by taking the limit of number m of periods going 
to infinity. We obtain from Eq. IT^ 



lim 

-n — >oo 



'At)) 



Ar(T) 



1 



-r(r) 



(15) 



It is important to notice that the steady state value de- 
pends on the duration r of the periods. Moreover, since 
A and 7 depend on r, consequently also F and Ap de- 
pend on r [c.f. Eqs. ^ and ([9])], and the steady state 
value may also depend on the form of the environmental 
spectral density applied in each period, and not only on 
its temperature T. 

As a cross check, it is easy to see that with increas- 
ing duration of the period r, we obtain correctly the re- 
sult without shuttering. The denominator in Eq. (jl5p . 
1/ (1 — e"'"*-'^-'), goes to unity with increasing r and we 
obtain {n)s = Ar(T) — 1/2. This matches with the long 
time result without the shuttering which can be calcu- 
lated by taking the limit t — > c» in Eq. ([7]). 

To gain more physical insight, we expand exp [— r(r)] 
which appears explicitly and implicitly in Eq. psp . with 
respect to F. Keeping the terms to lowest order in F, and 
neglecting the term 1/2 since we consider high tempera- 
ture case, we obtain 



{n)s 



l Iodt'Ajt') 
2f,dt'^{t')- 



(16) 



Moreover, this equation is valid also in the Markovian 
limit since with the Markovian values for A and 7 we 
obtain (n)^ = k-^T /luq, which matches the high temper- 
ature Markovian result [33 . 
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Equation ((16)) demonstrates that the ratio between the 
time averaged diffusion and dissipation gives the steady 
state value of the average system energy in the presence 
of shuttered reservoirs. Here, the average is taken over a 
single period of duration r. The system is not in ther- 
mal equilibrium with its environment and has an effective 
temperature Tcff which is different than the temperature 
T of the environment. Moreover, the steady state average 
energy of the system can be controlled by the period du- 
ration T and the environment " on-resonance" parameter 
r. 

The steady state of the system is a thermal state, see 
Eq. (fT3|) . with (n(i)} = (n)s- The corresponding effective 
temperature can be written as 



Toff = !—("-) 6 

Kb 



(17) 



This differs from the environment temperature T since 
{n)s 7^ n. Moreover, since the system reaches an ef- 
fective thermal state, it eventually fulfills the detailed 
balance condition with appropriate modifications to the 
transition rates given by the shuttered reservoir. 

We have presented above the formalism and the cross- 
checks for its validity. In the following section we con- 
tinue with the results for the heating function dynamics 
from Eq. p4|) and the steady state values of the average 
system energy from Eq. (|15p . 



IV. THE SYSTEM DYNAMICS 
A. Short and intermediate time dynamics 

Figure [T] displays the heating function dynamics over 
a couple of first shuttering periods. It illustrates how 
the system exhibits (a) Zeno and (b) anti-Zeno effect. 
For various types of non-Markovian dynamics with tra- 
ditional unshuttered reservoirs in this regime see Fig. 1 in 
Ref. [i| and for a study of the Zeno - anti-Zeno crossover 
see Ref. Q- The key aspect for this paper is that both 
effects are possible for short evolution times. Moreover, 
one can control which one appears by controlling the pe- 
riodicity T and the properties of the environment. Note 
that the figure includes also the non-Markovian dynam- 
ics within each period and not only the coarse graining 
in T [c.f. Eq. dm)]. 

Examples of the heating function dynamics for inter- 
mediate times, when the system approaches the steady 
state, are displayed in Fig. O The results demonstrate 
clearly that dynamics which begins with reduced heating, 
and therefore as Zeno- type for short evolution times, be- 
comes anti-Zeno for long times. This can be observed as 
a crossing between the dotted and solid lines in Fig. [5] (a) 
which shows in detail the passage from Zeno to anti-Zeno 
region. 

Moreover, Fig. [21(b) demonstrates that the approached 
steady state value for the dotted line is around 3.5 times 
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FIG. 1: (Color online) The short time behavior of the heating 
function displaying (a) Zeno and (b) anti-Zeno effect. In the 
Zeno effect, the shuttered reservoir reduces the heating of the 
oscillator while for the anti-Zeno case the heating is enhanced 
when compared to the dynamics with the traditional unshut- 
tered reservoir. The parameters are g — 0.1, r — 0.5/u>c, 
n = ksT/ujo = 10. In (a) r = 10 and (b) r = 0.1. 



the one expected from the temperature of the environ- 
ment. Furthermore, it is striking that there exists a cross- 
ing between the dotted and dashed lines in Fig. [5] (b) . In 
another words, the dynamics which starts as Zeno- type, 
shows for long times more pronounced anti-Zeno charac- 
ter than the case which already initially begins as anti- 
Zeno type. As a consequence, the highest steady state 
values of the energy for long times are reached for the 
parameters that in short times give QZE and reduced 
decay. 

The results above clearly demonstrate that the appear- 
ance of QZE and AZE is time-dependent in QBM. There 
exists a crossover between the two effects in time. The 
fact that the appearance of these effects depends on the 
total evolution time of the system, in addition of the pe- 
riodicity T, is in contrast to more simple systems like 
two- level atoms. It has been shown in Ref. 0], which 
studied the controlling of the Zeno and anti-Zeno ef- 
fects in QBM for short times, that the appearance of 
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FIG. 2: (Color online) The heating function dynamics for in- 
termediate times when the system appproaches a steady state. 
In panel (a) ^ cOct ^ 0.15 x 10^ which shows the small time 
blow up of panel (b) where ^ ujct ^ 10'*. The parameters 
correspond to those used in Fig.[T]except r — 1/ujc- The solid 
line is for an unshuttered reservoir showing the thermaliza- 
tion to {n)/n = 1 in panel (b). The dashed line is for the 
case of pure anti-Zeno dynamics by shuttered reservoir for all 
times. The dotted line is for shuttered reservoir when the 
dynamics begins as Zeno-type (reduced heating) but crosses 
later to anti-Zeno behavior (increased heating). This crossing 
between the solid and dotted lines, which occurs at uJct ~ 800, 
is shown in detail in panel (a). 



the QZE and AZE depends on the quantum number n 
for QBM. Even though the time averages of the effec- 
tive decay coefficients remain fixed, this n-dependence 
allows the crossover from Zeno to anti-Zeno effect to oc- 
cur with time. When the system heats, the expectation 
value (n) increases with time, and the system experi- 
ences a crossover from Zeno to anti-Zeno region. This 
is not possible for a two-level system and demonstrates 
that the Zeno~anti-Zeno phenomena has more variety in 
QBM than in more simple systems. 



B. Long-time dynamics and steady state 

The results presented in the previous subsection sug- 
gest that the heating function with shuttered reservoirs 
approaches asymptotically values which are higher than 
with the unshuttered reservoirs. Actually, this is also in- 
dicated by the analytic solution given in Eq. I|14p where 
the denominator is always between and 1. The sys- 
tem which is forced repeatedly to go through the non- 
Markovian behavior, experiences enhanced dissipation 
compared to the Markovian one. The system reaches 
a steady state since the time dependent diffusion and 
dissipation coefficients, A and 7 respectively, have well 
defined time averaged values over t. 

By using Eq. ([TS]) one can calculate the steady state 
system energy as a function of r, r, and the coupling con- 
stant g. The magnitude of the coupling constant, in the 
region of the validity of the model (weak coupling g <ti I) 
does not affect the steady state energy. Moreover, {n)s, 
as a function of r, varies in a rather small range for r < 1 
and then reaches a stable value for increasing r. Thus, 
the most interesting parameter here is the shuttering in- 
terval T. This conclusion can be reached also by making 
a series expansion of {n)s, with respect to small t. The 
dominant term contains only r, and the parameters r 
and g come into play only in the higher order terms. 

Figure [3] displays the steady state (n)s as a function of 
T, where for other parameters some convenient values has 
been chosen. As expected, for increasing r, the average 
energy of the steady state approaches the one with un- 
shuttered reservoir and corresponds to the temperature 
of the environment. When the duration of the periods 
T is decreased, the steady state energy of the system in- 
creases rapidly. This can be explained by studying the 
behavior of A and 7 as a function of time. For the typ- 
ical spectral density we use in this paper [c.f. Eq. (|4])], 
A initially increases in the non-Markovian region faster 
than 7 [c.f. Eqs. ^ and The feature is enhanced 

in the time averaged values of these quantities when r 
decreases and demonstrates quantitatively the behavior 
seen in Fig. [31 

Naturally, the energy conservation has to be satisfied, 
ft requires work to periodically shutter the reservoir. In 
other words, it takes work to periodically couple the sys- 
tem with the reservoir. This can be seen as the origin of 
the increased system energy and the energy is conserved 
since the minimum work required exceeds the oscillator 
energy (soj . 

V. DISCUSSION AND CONCLUSIONS 

The creation of artificial environments and the reser- 
voir engineering techniques open new avenues to study 
the non-Markovian dynamics of open quantum systems 
and to improve the capability for the control of quan- 
tum systems. We have taken here an initial step towards 
time-dependent reservoir engineering by considering the 
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FIG. 3: The steady state value of the heating function (n)s as 
a function of the shuttering period r. Parameters are a = 0.1, 
r = W,n = 10. 

case when a structured reservoir is switched on and off 
in a periodic manner (shuttered reservoir). 

The results for quantum Brownian motion show that 
the interaction of the reduced system with the shuttered 
reservoir changes the system dynamics drastically com- 
pared to the conventional dynamics. For short times the 
system may exhibit Zeno or anti-Zeno behavior depend- 
ing on the properties of the environment. A striking 
dynamical feature arises for long-times - the system al- 
ways reaches a steady state in which the average system 
energy is larger than the one corresponding to the the 
temperature of the environment. The shuttered reser- 
voir forces the system to experience the non-Markovian 
behavior periodically and the steady state properties are 



consequently given by the time averaged properties of 
the diffusion and dissipation between successive reservoir 
shuttering events. 

One of the interesting consequences is that the appear- 
ance of Zeno or anti-Zeno effects becomes dependent on 
the total evolution time of the system. What starts as a 
Zeno effect for short and intermediate times turns to anti- 
Zeno effect when the steady state regime is approached. 
The Zeno studies in simple systems, such as a two- level 
system, show how either Zeno or anti-Zeno effect appears. 
Here, we have demonstrated that a crossing between the 
two effects occurs with time and in QBM the Zeno phe- 
nomena is very rich. 

In the future, advanced reservoir engineering tech- 
niques may allow new ways to implement quantum con- 
trol. We have shown here that shuttered reservoir gives 
rise to interesting dynamical effects. A next step is to 
study the case where the properties of the environment 
remain the same for a given "on-period" but vary be- 
tween different periods. Very precise ways to control 
the system properties by using reservoir engineering com- 
bined with laser cooling may open new possibilities for 
the development of general quantum simulators for open 
systems whose dynamics is traditionally tedious to solve 
by analytical means. 
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